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ABSTRACT 


Biaxial  am)  uniaxial  cxper iments  have  been  conducted  on 
a  thin  sheet  of  natural  rubber,  which  can  be  assumed  to  be  incompress¬ 
ible,  isotrcpic,  and  perfectly  elastic.  The  strain  energy  function  and 
constitutive  equations  have  been  determined,  and  the  material  is  c'las- 
si  ied  as  a  Generalized  R i v? i n- Mooney  type. 

B  la;:  I  a  I  experiments  were  then  conducted  on  the  same  sheet  with 
a  circular  cutout  and  stress  concentration  factors  were  obtained. 

Results  indicate  a  significant  increase  in  the  factor  with  increased 
d i splacements . 

A  modified  Particle  -  in  -  Cell  (P.l.“.)  method  has  been 
developed  and  analytical  results  ware  obtained  for  a  sheet  with  a 
rigid  circular  inclusion,  it  is  shown  that  the  stress  concentrat ion 
factor  for  a  Rivlin-Mooney  material  increases  with  increasing  deforma¬ 
tions,  a  result  which  is  in  qualitative  agreement  with  solutions  ob¬ 
tained  by  other  methods  [31»32j.  The  use  of  the  Generalized  Rivlin- 
Kooney  material,  however,  leads  to  a  decrease  in  stress  concentration 
with  increasing  deformations. 


LIST  OF  SYMBOLS 


contravar lent  components  of  in-plane  metric  tensors 
of  undeformed  and  deformed  body 

Rivl in-Mooney  material  constants 

contravar ient  components  of  metric  tenors  of  undeformed 
and  deformed  body 

strain  invariants 

stress  concentration  factor  at  boundary  of  rigid 
inclusion  (ratio  of  radial  stress  at  boundary  of 
inclusion  to  that  at  edge  of  sheet) 

stress  concentration  factor  at  hole  boundary  (ratio  of 
tangential  stress  at  edge  of  hole  to  that  at  edge 
of  sheet) 

cell  spacing  used  in  °.I.C.  method 
stress  resultant 
shear  resultant 
radius  of  undeformed  hole 
displacement  components 
displacements  at  edge  of  sheet 

strain  energy  density  function  (per  unit  volume  of 
unde^ormed  body) 

coordinates  of  a  point  in  deformed  body 
coordinates  of  sides  of  deformed  sheet 
coordinates  of  a  point  in  undeformed  body 


coord! nates  of  sides  of  tndeformed  sheet 


cr  yi 


'j 
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X 

y 


(or  >  ,), 
(or  \g) 


physical  component  of  stress  tensor 

ratio  of  deformed  to  undeformed  length  of  membrane 

elongation  ratios  in  x  (or  1)  and  y  (or  2)  directions 
respectively 


('•)„»(*  )w  ratio  of  deformed  to  undeformed  lengths  in  x  and  y 
y  directions  respectively 

~  J  component  of  stress  tensor 


'  'fc 


I.  INTRODUCTION 


An  ever  increasing  use  of  reinforced  and  unreinforced  poly¬ 
meric  materials  in  complex  technological  situations  has  given  impetus 

to  the  general  study  of  elastomeric  materials  subject  to  large  deforma- 
r 


t  ions 


But,  owing  to  the  nonlinearities  of  the  basic  equations 


for  such  materials,  analytical  solutions  have  been  obtained  for  only  a 
limited  number  of  problems.  in  general  the  material  is  assumed  to 
be  perfectly  elastic  and  incompressible,  and  the  deformations  are 
homogeneous  or  highly  symmetric,  so  that  nonlinear  partial  differential 
equations  can  be  avoided  in  the  analysis  (see,  for  example.  Chapter  III 
of  Ref.  7  for  discussion  of  the  work  of  Rivlin,  also  see  Refs.  6,  11-18), 

Although  the  literature  dealing  with  the  classical  stress 
concentration  problem  for  linear  elastic  materials  abounds,  (see  for  •**« 
example  Refs.  19-22)  only  very  few  investigations  can  be  found  which 
realistically  treat  the  stress  concentration  of  elastomeric  materials 
subject  to  large  deformations.  The  reason  lies  in  the  difficulty  in 
solving  the  nonlinear  partial  differential  equations  which  govern  the 
stress  and  displacement  fields. 

One  way  of  attacking  such  problems  is  to  use  the  method  of 
successive  approximation  to  obtain  solutions  In  such 


cases  the  stress  and  displacement  fields  are  expanded  into  power  series 
in  a  characteristic  real  param&ter  e,  which  depends  on  the  nature  of 


the  problem.  When  these  series  expansions  are  substituted  into  the  field 
equations  of  finite  elasticity,  the  coefficients  of  the  successive  powers 
of  s  in  the  resulting  identities  yield  linear  equations  for  the  deter¬ 
mination  of  the  successive  terms  in  the  power  series  expansions.  The 
parameter  £  is  so  chosen  that  the  identities  corresponding  to  the  first 
degree  in  e  coincide  with  the  set  of  field  equations  in  the  classical 
linear  theory.  Because  of  the  lengthy  calculations  required  most  of  the 
problems  actually  solved  in  this  manner  have  been  limited  to  the  first 
and  the  second  degree  in  s.  In  particular,  Adkirtet  alL  applied 

this  method  to  plane  problems,  and,  by  using  the  complex  variable  technique 
obtained  second-order  solutions  for  an  infinite  plate  with  a  circular 
hole  or  a  circular  rigid  inclusion,  subjected  to  a  uniaxial  tension  at 
infinity.  The  results  showed  that  for  a  Mooney  material  the  degree  of 
stress  concentration  at  the  hole  is  reduced  as  the  load  is  increased. 
Similar  observations  were  made  by  Guz1,  Savin  and  Tsurpal1  in  the<r 
study  of  an  infinite  plate  of  physically  nonlinear  material  having  a 
circular  hole.  They  assumed  that  the  displacement  field  was  within  the 
range  of  applicability  of  infinitesimal  strain  theory.  In  a  recent  paper, 
Evans  and  Pister^^  investigated  the  same  problem  by  using  a  power  series 
expansion  up  to  the  third  order.  They  showed  that  unless  the  deviation 
from  linearity  in  the  stress-strain  relation  is  very  small,  the  series 
solutions  do  not  converge  rapidly  enough  to  lead  to  reasonable  approxima¬ 
tions  to  the  exact  solutions. 


3. 


Another  way  of  solving  the  resulting  nonlinear  equations  of  the 

stress  concentration  problem  is  to  employ  numerical  procedures.  Using  a 

forward  integration  method,  Rivlin  and  Thomas^^  solved  the  problem  of 

the  axially  symmetric  deformation  produced  in  a  circular  sheet  of  a  Mooney 

material,  containing  a  central  circular  hole,  when  subjected  to  radial 

tensile  forces  at  its  outer  periphery.  Their  displacement  results  were 

found  to  be  in  good  agreement  with  the  experimental  measurements  on  a 

test  piece  of  vulcanized  natural  rubber.  No  results  on  stress  concentra- 

[321 

tion  was  presented.  Recently,  Yang  re- invest igated  the  axisymmetric 
problem  and  was  able  to  reduce  the  equations  to  two  decoupled  first-order 
equations,  which  were  solved  numerically  by  the  Runge-Kutta  method.  The 
results,  for  both  the  hole  and  the  rigid  inclusion,  indicate  that  the  stress 
and  strain  concentration  factors  increase  with  increasing  load  for  the 
Mooney  material.  An  analysis  of  the  displacements  obtained  by  Rivlin  and 

f3 1] 

Thomas'  reveals  the  same  trend  for  the  concentration  factors. 

n 

In  this  thesis,  the  problem  of  a  highly  elastic  sheet  with  a 
centered  circular  hole  or  rigid  circular  inclusion,  subjected  to  a  general 
biaxial  stress  state,  is  considered.  Since  we  consider  very  large  deforma¬ 
tions  and  general  biaxial  stress  states,  none  of  the  previous  approaches 
are  applicable.  Instead  the  nonlinear  partial  differential  equations  are 
linearized  by  adopting  an  incremental  procedure.  A  variation  of  the 
Part icle- In-Cel  1  (P.I.C.)  method,  which  has  been  so  successfully  applied 
to  hydrodynamic  blast  problems1  ,  is  developed  specifically  for 


k. 


this  type  of  elliptic  problem.  Essentially,  the  method  combines  the 
features  of  the  Euler ian  and  Lagrangian  approaches  by  introducing  an 
tuleri:.n  computational  mesh  and  by  keeping  track  of  the  material  particles 
as  in  a  Lagrangian  formulation.  A  description  of  this  method  can  be 
found  in  Chapter  II,  Section  H. 

Two  series  of  experiments  have  been  conducted.  One,  was  carried 
out  on  the  sheet  without  the  hole,  in  order  to  characterize  the  mechanical 
properties  of  the  rubber  for  use  in  the  analytical  calculations.  The 
second  set  were  the  biaxial  stress  concentration  experiments  performed  on 


the  sheet  with  the  hole. 


5. 


2.  THEORY 


A.  Equilibrium  Equations 


for  a  Thin  Plane  Sheet 


[7,10] 


Consider  a  flat  elastic  sheet  whose  thickness  before  deforma¬ 
tion  h  is  constant  and  is  small  compared  to  its  other  dimensions, 
o 

The  sheet  undergoes  a  finite  deformation  symmetrical  about  the  middle 
plane  z  =  0  which  thus  becomes  the  middle  plane  Z=0  in  the  deformed 
state.  X,  Y,  Z  are  the  coordinates  of  a  point  in  the  deformed  body 
which  before  deformation  was  at  x,y,z,  where  both  sets  of  coordinates 
are  referred  to  the  same  rectangular  cartesian  system. 

The  equilibrium  equations,  in  the  absence  of  body  and  inertia 
forces,  are: 


_  'J 


(2-1) 


w’lere  the  double  line  denotes  covariant  differentiation  with  respect  to 
i(X,Y,  or  Z).  Since  in  rectangular  coordinates  the  components  of  the 
stress  tensor  r 1 coincide  with  the  physical  components  of  the  stress 
c..f  Eq,  (2-1)  can  be  written  as: 


a 


'J  ' 


=  0 


(2-2) 


where  a.. 

'J 


is  the  stress  per  unit  area  of  the  deformed  body. 


6 


Integrating  the  equilibrium  equations  (2-2),  and  noting  that 
the  surfaces  of  the  plate  are  traction  free,  we  obtain 


W°  fo's  =  XiY) 


(2-3) 


where 


"of  ■  / 


-h 


*06  dZ 


(2-4) 


is  the  stress  resultant  per  unit  length  of  the  deformed  plete  and  h  is 
Its  thickness.  The  shear  resultant 


q^_  =  f  r  dZ  =  0 
Ma  Gts. 


(2-5) 


due  to  the  symmetrical  stress  distribution  about  the  middle  plane. 

We  define  a  scalar  function  X  which  represents  the  extension 
ratio  at  the  middle  plane  Z=0  in  a  direction  normal  tc  it, 


Thus 


<z'°> 


h 


o 


=  /,  im 
Z-*  0 


(2-6) 


(2-7) 


From  (2-4): 


n  .  ,  h 

cae  i  r  j7 

2h  =  2h  J  aoe  dZ 

o  o  -h 


(2-8) 


and  from  (2-7) 


h  -♦  0  o 
o 


Since  the  faces  of  the  sheet  are  free  from  applied  forces 


/rsf-zz’z.o-  ° 

o 


For  very  small  hQ  we  can  approximately  write 


t°zz]z  =  0  *  0 


"06  '  2ho  ^(<Wz=0 


(2-9) 


(2-10) 


Substituting  (2-10)  into  (2-3),  we  obtain  the  final  form  of  the  membrane 
equilibrium  equations 


0  <2*°» 


(2-11) 


5x  (Hx)  +  57  f'Xc!xr  =  0  (z  *  o> 

57  (XaXY}  +  57  (XaYY)  -  0  <Z  *  0) 


(2-12) 


B.  Constitutive  Equations 


If  the  sheet  Is  made  of  a  homogeneous,  isotropic,  elastic  and 


Incompressible  material 


V  1 


(2-13) 


and  the  strain  energy  density  function  W  (measured  per  unit  volume  of 
the  undeformed  body)  takes  the  form 


W  =  w(Tr  t2) 


(2-14) 


The  inplane  stresses  on  the  middle  plane  can  then  be  written  as^0^ 


ro3 =  Cce  3  H  aQS  +  L  aQS  (a,a  =  X»Y) 


(2-l*o) 


where 


”  =  2lrr  +  >2$;> 


L  =  -  2\ 


2(\2T1  -  1A2)  M- 


(2-16) 


.  Q6  .00  _  .  , 

and  a  *  A  are  the  'np  ane  components  at  Z =  0  of  the  metric  tensors 
(g1-^)  of  the  undeformed  and  (GIJ’)  deformed  body,  respectively.  Thus 
the  general  metric  tensors  at  Z = 0  are  written  as; 


9. 


10. 


and 


(2-21) 


where  >  -  X(X,Y)  represents  the  extension  ratio  at  the  middle  plane 
in  the  Z  direction  and  is  assumed  to  be  finite  and  nonzero.  We  now 
can  substitute  Eqs .  (2-20)  and  (2-21)  into  (2-17)  and  (2-18)  to  obtain 
the  values  of  the  components  of  the  metric  tensors  g ^ at  the 
plane  Z  =  0. 


9XX  “ 

axx  ■ 

(*1>2 

5xy  = 

aXY  = 

3x 

3X  hY 

SX  7>Y 

3yy  * 

aYY  = 

(~)2  + 

9ZZ  = 

(If)2 

CSC 

X 

M 

It 

9yz  x 

e 

or  g, . 

U 


aXX  *XY 


aXY  aYY 


0  l/\c 


(2-23) 


and 


g  =  Ig.jl  =  aA' 


(2-24) 


where 


a  =  ^  ao8 ■  "  aXX  aYY  '  ^aXY^ 


(2-25) 


The  component  of  the  contravar iant  metric  tensors  g'j  and  a08  can 
be  determined  from  the  equation: 

j.  cofactor  g{| 

g  =  - - - L 


(2-26) 


and  are  therefore 


XX 

x> 

g 

=  a  = 

XY 

XY 

g 

=  a  = 

YY 

YY 

9 

=  a  a 

22 

a 

9 

a  A2 

X2 

Y2 

<J 

=  g  =  i 

XX 

XY 

a 

a 

XY 

YY 

9 

a 

0 

0 

Ig'-M 

=  \2  [a 

axyA" 

2 


,2  ~  a 


(2-27) 


(2-28) 


(2-29) 


12. 


The  strain  invariants  art  given  by  the  expressions 


,  rs  _  .  2  rs  . 

T.  =  g  G  =  a.  +  a  A 

rs  rs 


? 


r  „rs  t  1  . r s 

T2'9rsG  T3  =  ^+arsA 

.2, 

T3  =  G/9  =  ^=  . 


(2-30) 


and  therefore  a  =  \  because  of  incompressibility, 

2 

aw  +  a\/3/  +  * 

Finally,  then 


t  =  aXX  +  aVY  +  3 


1 


T2  =  SXX  +  SYY  +  Wa 


(2-31) 


The  most  general  power  series  form  of  W  for  an  incompressible 

r  T7i 

homogeneous  Isotropic  elastic  material  is 

00  00  . 

w  =  v:  E  c..  (T  -3)1  (T_-3)J  (2-32) 

i=0  j=0  J  1 

where  C..  are  constants  and  C  =0. 
i  j  oo 

In  the  following,  we  assume  a  R i vl i n-Mooney  type  of  material 
which  is  expressed  by 

W  =  C1(T1-3)+C2(T2-3)  (2-33) 

The  necessary  and  sufficient  conditions  for  W  >  0  is  that  >  0  and 
Cp  >  0.  Substitution  of  (2-31)  and  (2-33)  Into  (2-16)  yields 


i 


(2-34) 


L  *  —  2a  Cj  —  2  (a^  +  ^YY  ~  ^^^2 
Finally,  the  stresses  of  Eq.  (2-15)  can  be  written  as  follows: 

°XX  “  2CI  ~t  ■  2C2  aXX  '  2aCl  +  2  ~a 

avx  C? 

°rv  -  2CI  -r  *  2C2  aYY  •  2aCl  +  2  T  (2-J5) 

°XY  *  •  2(CI  +  aC2)  "T 

C.  Displacement  Equilibrium  Equations 

Replacing  X  and  Y  by  the  indices  1  and  2,  respectively, 
and  appealing  to  Eqs.  (2-35) »  we  have 

\axx  =  *0| ,  =  2  Ta*  a22  "  2C2v/®’  al  1  “  2'^Cl  +  2 

C  a  C 

\ayy  =  \a22  =  2  -  2C2  a22  -  2\^C,  +  2  (2-36) 

a,2 

XaXY  =  Xa12  =  "  2(C1  +  aC2)  Va1 

The  equilibrium  equations  (2-12)  then  may  be  expressed  as 

2 

C[2a(aJ2,2  -  a22>J)  +  (a22  +  3a  )  a^  -  aJ2  a,2] 

2 

+  a, |  +  aCa^a.j  +  a]2  a>2)  +  2a  (a,2,2  +  a,,,,)  «  0  (2-37) 


I 


2 

snd  £f23(3|2*j  *  3j  j  ig)  +  (3j  |  3a  )  ^»2  *  a]2  a,j^ 

(2-38) 

+  a tg  *  a(a22  a.2  +  a|2  a,,)  *  2a2(a|2,,  a-  a^..,)  =  0 

where  C  =  Cj/C2  ,  (2-39) 

and  a  comma  denotes  partial  differentiation  with  respect  to  1  or  2 
(X  or  Y). 

We  define  the  displacements  in  terms  of  the  coordinates  of  the 
deformed  body 

u(X,Y)  =  X  -  x  (X , Y)  ‘ ' 

v(X,Y)  =  Y  -  y (X,Y)  (2-40) 

Thus  in  terms  of  these  displacements 

3jj  =  (1-u.j)2  +  (v,j)2 

a22  =  (u«2^2  +  (*-v»2^2  (2-41) 

a,2  =  "  [0-u,,)u,2+  (1-v,  )v(1] 
a  =  [ (l-u,j) (l-v,2)  -  u>2  v,}]2 


and  the  equilibrium  equations  in  terms  of  the  displacements  are 


t 


•  5- 

2  2 

(Vj^Vij  )  (u »22 * 2+v »22u *2  ^  ^  |  ^ *2"^2  1  ^  1 2Y2^*Y *  1 2^*2^ 

2  2  I) 

+  (y  +u>2)(u,11v2+v,11u,2)  +  3(y1Y2-u,2v,1)  (u,1jy2+vH2y1+u,21V,1+U,2v,11^ 

2  2  2 

+  ^u*iiV2+v*21V1+U,21  V,1+U’2V,lP  "  ^V1V2*"U'2V*  1 J  ^u,11^v*1  V2_3Y1V2+2u,2V’  1V1  * 

2  3  2  2  2  2 

+  U»12(-V,V’1“V*|  f2u'2V1V2+V,1V2_U,2  V,P+  u’22^2u’2V,1VrYlV2+v,l  V2^ 

2  2  3  2  2  2 

+  v*  j  j  (“V ^  u,2^*2 v 1 1  V] V2“3 u »2 v *  j  )+  v » |2  (“v i ”v*i  v  1+0*2  V1+2u’2'/,1V2“V1V2^ 

+  v*22(u,2V2v’lVr2“U,2V,l2^  =  °  (2-42a) 

and 

2  2 

C[  (v 2+u>2  )(v>  j  | ^ j  j  v»  j )  +2  (%2v » |  |  u »g)  (^»  |2^ 1 2^ ^ 

2  2  4 

+  (vj+v»j  )  ( v,22^1  1  ^  ^  (^2^|" i  ^ >2^  ^>22"^2^>  12^*^>  12^,2^*^,  1  ^*22^ 

+  (v»22v i+u*  12V2+V*  12U,2+V’  1  U,22)  ~  (V2''l"V' 1  U’2)2[ V’22(*U^Vr3YlY2+2v’ 1  U*2V2) 

2  3  2  2  2  2 
+  v»i2^“U,2V2”U,2  f2v’lVlV2+U,2VrU,2V,l  v*  j  ^2u,?v*1V2"V1V2+U,2  Yr 

2  2  3  2  2  2 
■*■  u >22 ("v2 v*  1 +2u ’2 ^  1  V2~3 11  ’2  ^ * i  2 'Y2"**^,2  Y*2"^"'^*l  "''2"*"2u,2v,1Y1_Y1Y2^ 

2  2 

+  u»1 1 (v,jV2+2u,2YjY2-u,2  v, 1 ) ]  =  0  (2-42 b) 

where  Vj  =  1  -  u,j  ,  Yg  *  1  -  v>2  (2-43) 

D.  Piecewise  Linearization  of  the  Equilibrium  Equations 

The  two  equilibrium  equations  (2-42)  are  second  order  nonlinear 
partial  differential  equations  in  u  and  v.  In  order  to  solve  for  u 
and  v  we  assume  that  the  body  is  strained  to  its  final  position  by 


successive  steps  of  small  displacement  increments  Au  and  Av.  The 
total  displacement  after  n+1  such  steps  is  expressed  as 
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(2-44) 


where 


n 

n  „  ,  i 

u  =  l  Au  ,  v 

i=l 


n 

r 

i=l 


(2-45) 


It  is  assumed  that  Au,  Av  are  sufficiently  small  so  that 
terms  containing  squares  or  higher  order  terms  in  Au,  Av  or  their 
partial  derivatives  can  be  neglected  with  respect  to  first  order  terms. 
Thus,  for  example, 


/,  n+1 .2 

(1-u,,  ) 


„  n  ,  n+1  2 
(l-u.j-Au.j  )  1 


/,  n  .2 

( 1  -  u ,  ^ ) 


-  2(l-u?1)Au?|1  +  0(; 


n+1 

u,, 


(2-46) 


Following  this  assumption  all  the  expressions  appearing  in  the 
equilibrium  equations  are  rewritten  by  expressing  u,  v  in  the  forms  of 
(2-44)  and  retaining  only  first  order  terms  in  Au,  Av  and  their  deriva¬ 
tives.  For  simplicity,  the  superscripts  of  u11,  v11,  Aun+^ ,  Av1^  will 
be  omitted  in  writing  the  following  linearized  form  of  the  equilibrium 
equations  at  the  (n+l)^  step 


1 


^11  ^  ^22  Au*22  ^12  ^u*  12  ^2  ^*2 

+  B j  ^  iv,  j  j  +  ®22  Av»22  ®J2  A^*|2  "*"  ^1  A^*]  ®2  A^*2  ^|  “  ^  (2*47s) 

and 

Dj  |  1 1  +  t>22  ^u,22  ^12  ^u*  12  *  ^1  ^^*i  ^2 

+  E^  ^v»n  +  E22  ,iv,22  +  E12  Av,12  +  E1  Av*l  +  E2  Av*2  +  *2  ~  0  (2-47b) 

The  coefficients  A^,  ....  . Kj ,  . j . K  are  to 

be  found  in  Appendix  A  and  are  functions  of  the  derivatives  of  the  total 
displacements  at  the  completion  of  the  n*^1  step  and  contain  the  RIvl  in- 
Mooney  constants. 

E.  Modified  P.I.C.  Method 

The  linearized  equilibrium  equation?  (2-47)  are  solved  numeri¬ 
cally  by  using  a  modified  form  of  the  P.I.C.  Method,  which  has  been 

.  [34] 

successfully  used  in  wave  propagation  problems  .  The  procedure 
fol 1 ows . 

First  a  fixed  grid  (Eulerian)  is  introduced,  dividing  the  region 
into  a  finite  number  of  cells.  The  displacement  of  each  cell  is 
represented  by  the  displacement  at  its  center.  The  partial  derivatives 
of  Au,  Av  are  expressed  in  finite  difference  form  at  each  of  the  cell 
centers,  which  act  as  the  pivotal  points.  The  substitution  of  these 
finite  difference  forms  into  the  differential  equations  (2-47),  produces 
a  set  of  linear  algebraic  equations  in  the  set  of  unknowns,  Au  and  Av, 
at  each  point.  The  incremental  displacements  determined  in  this  first 
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(Eulerian)  part  satisfy  the  equilibrium  equations  in  the  sense  of  class¬ 
ical  linear  theory,  where  ;t  is  assumed  that  the  final  coordinates  are 
identical  with  the  initial  ones.  The  geometric  nonlinear  effect,  due  to 
changes  in  the  configuration  of  the  material,  is  accounted  for  in  the 
second  part  of  the  procedure. 

We  begin  by  distributing  particles  in  each  cell,  where  a  cell 
(i,j)  is  the  square  around  mesh  point  (i,j)  bounded  by  the  lines  i ±  k/2 
and  jik/2.  These  particles  move  during  each  step,  and  their  incremental 
displacements  are  assumed  to  be  the  same  as  that  of  the  cell  center  in 
which  they  find  themselves  at  the  beginning  of  a  step.  The  position 
and  total  displacements  of  these  particles  are  recorded  throughout  the 
deformation.  In  solving  for  the  incremental  displacements  of  the  n+l**1 
step,  for  example,  it  is  necessary  to  know  at  each  cell  center  the  total 
displacements  accumulated  during  the  first  n  steps.  These  displace¬ 
ments  ere  assumed  to  be  equal  to  the  average  displacements  of  all  the 
particles  in  the  cell  at  the  end  of  the  nth  step.  If  no  particles 
move  in  or  out  of  a  cell  then  its  displacement  remains  unmodified. 

At  the  end  of  each  load  step  the  nonlinear  equations  are 
checked  so  as  to  be  sure  that  rhe  assumption  of  piecewise  linearity  is 
reasonable.  If  the  check  fails,  the  step  has  to  be  refined  and  the 
computation  repeated. 

F.  Extension  of  a  Thin  Elastic  Sheet  with  a  Hole  at  its  Center 

Consider  a  rectangular  sheet  subjected  to  biaxial  extensions 
along  its  sides.  The  sheet  contains  a  traction  free  hole  at  its  center. 


The  method  of  solution  described  can  be  applied  to  a  hole  of  arbitrary 
shape.  If  the  hole  is  symmetr ical  with  respect  to  the  axes  of  symmetry 
of  the  sheet,  then  it  is  sufficient  to  consider  only  one  quarter  of  the 


sheet.  in  what  follows,  a  hole  with  such  synretry  is  assumed,  and  the 
analysis  is  thus  carried  out  for  only  one  quarter  of  the  sheet.  The 
surface  of  the  sheet  is  divided  into  cells  by  introducing  a  grid  system 

(Fig.  i).  The  pivotal  points  are  located  at  the  cell  centers  and  have 

a  constant  spacing  k.  The  partial  derivatives  at  i , j ,  expressed  in 
central  difference  form  are: 

u,|  (i,j;  =  [u(!+k,J)  -  u(i-k,j)]/2k  +  0(k2) 

=  [u(i,j+k)  -  u(  i  ,j-k)]/2k  +  0(k2}  (2-48) 

u,n(i,j)  =  [u(i+k,j)  -  2  u  ( i ,  j )  +  u(i-k,j)]/k2  +  0(k2) 

u,22(i,j)  =  ( j(i ,  *+k)  -  2u(i ,  j )  +  u(i,j-k)]/k2  +  0(k2) 

’J  ^  =  i  u(i+k,j+k)  -  u(i-k,  j+k)  -  u(i+k,j-k)  +  u(i-k,j-k)] /4k2  +  0(k2) 

introducing  ^uch  forms  for  the  partial  derivatives  into  the 
equilibrium  equations  (2-47)  results  in  a  set  of  algebraic  equations 
with  each  equation  containing  the  unknowns  Au,  Av  of  point  (i,j)  and 
the  eight  surrounding  points.  Thus  the  equilibrium  equations 

in  central  difference  form  are 

A12  Au(i-k.j-k)  +  (4A22-2kA2)  Au(i .  j-k)-A12Au(i+k,j-k)  +  (4A]  ]-2kA1  )Au(i-k,j) 
-  8  (A  j .  +A22 )  Au(i,j)+ (4A]J+2kA1)  Au(  i+k,  j)  -  A]2iu(  i-k,  j+k) 

+  (4A22+2kA?)  Au(i,j+k)  +  A]2  Au(i+k,j-:k)  +  Z  Av(i-k.j-k) 


(continued  on  next  page) 
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+  (tS^-SkBg)  Av(i.j-k)  -  B12  Av(i+k,j-k)  +  (4Bn-2kB1)  Av(i-k,j) 

-  8<8H+S22)  Av(i'J)  +  (i+811+2kB1)  Av(i+k,j)  -  B]2  Av(i-k,j+k) 

+  (1»B22+2kB2)  Av(i,j+k)  +  BJ2  Av(l+k, j+’  +  Kj  =  0  (2-49a) 

and 

°12  +  (it022"2kD2^  &u( i ,  j-k)  -  0j2  Au(i+k,j-k) 

■*  (40n-2kD1)  Au(l-k.j)  *  8(Dn+022)  Au(i,j)  +  (4Dn+2kDJ)  Au( i+k, j) 

-  Dj2  Au(i-k,j+k)  +  (402O+2kD2)  Au(i,j+k)  +  0?2  Au(l+k,j+k)  +  E,2  Av(i-k.j-k) 
+  (^E22-2kE2)  Av(f.j-k)  -  E]2  Av(i+k,j-k)  +  (4EJ1-2kE1)  Av(i-k.j) 

-  3(En+E22}  A v ( i , j )  +  (4E]1+2kE1)  Av(i+k.j)  -  EJ2  Av(i-k,j+k) 

+  (4E22+2kE2)  Av(i,j-K':)  +  E12  V\v{l+k,j+k)  +  «  0  (2-49b) 

Boundary  Conditions 

The  x  -  and  y  -  axes  are  axes  of  symmetry,  and,  consequently, 
u(!,j)  *  0  on  the  y-axis, 

and  u(l-k,j)  =  -  u(l+k,j)  ;  v(l-k,  J)  *  v(l+k,j),  (2-50) 

while  v ( i , I )  *  0  on  the  x-axis, 

and  u(i ,  1-k)  =  u(  i ,  1+k)  ;  v(i,l-k)=-  v(i,l+k)  (2-51) 

Let  x  =  X|  and  y  =  y^  be  the  sides  of  the  undeformed  sheet  (Fig.  1) 
then  at  X  = 


21 


and 


u  =  (prescribed) 


n  =0 

XY 


(2-52) 


and  at  Y  =  Y, 


v  =  (prescribed) 


and 


'XY 


»  0 


(2-53) 


It  is  assumed  that  the  sides  X  =  X^  and  Y  =  Y^  are  far  enough  from 
the  hole  so  that  we  can  replace  the  shearing  stress  conditions  by  the 
linear  displacement  conditions, 


V  '  Y,  V1 


u  =  r  u  i 


(2-54) 


The  hole  is  free  of  traction  so  that 


XX 

’XY 


*  +  <j„„  m  =  0 


XY 

'yy 


l  +  m  =  0 


(2-55) 


where  the  direction  cosines  of  the  unit  normal  are 


v,  dX  dY 
x  =  cos (N,X)  -  dN  =  ds 


dY  -dX 
m=  cos(N.Y)  =  ~js 


(2-56) 


8y  substituting  Eqs .  (2-56)  into  (2-55)*  the  boundary  conditions  along 
the  hole  can  be  written  as 
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'XX  dS 


dX 
rXY  d‘J 


dY 
”XY  HS 


YY  dS 


(2-57) 


where  the  stresses  can  be  expressed  in  terms  of  the  derivatives  of  the 
displacements  by  utilizing  Eqs  .  (2 -?6)  and  (2-41).  The  piecewise  linear 
form  of  the  stress-displacement  relationship  can  then  be  expressed  as 


acxx=  a j  *. u ,  j  +  +  bl“V,l  +  b2^V,2  +  kl 


acYY  =  d,iu#l  +  d2*u,2  +  fijiv,,  +  e2iv,2  +  kg 
arrXY  =  fliu,l  +  f2Au,2  +  9li‘v,l  +  92^v,2  + 


(2  •  58) 


where  a j  $  a2 ■  b| t  b2 *  k ^  >  d ^  >  d^)  6| t  e2 >  kg  a nd  f| »  fg »  91*  92 *  ^ 
are  given  in  Appendix  A. 

dX  dv 

The  procedure  for  evaluating  ^  ,  ~  is  as  follows.  The 

hole  is  delineated  by  placing  evenly  spaced  particles  on  the  boundary 

curve.  These  particles  then  trace  the  boundary  throughout  the  deforma- 

d  X  d  Y 

tion  process.  The  derivatives  —  ,  —  at  a  particular  boundary  partic’e 

d  o  db 

i  is  obtained  by  expressing  X  =  X(S)  and  Y  =  Y (S )  at  the  points  i  -  l 
and  i  +  1  by  a  Taylor  series  expansion  about  point  i  (Fig.  2)  and 
retaining  only  terms  up  to  the  third  order.  The  solution  of  the  two 
simultaneous  equations  then  results  in  the  following  expressions 

dX,.,  -  X"-')<Si.!+l>2  +  ’‘(i>KSi,i+l)2  -  <r.-l,I)2J  .  .,.2, 

■ - h-i.,  (  1 


dY,.,  v(l+l)(s|_|  ,)  -  +  Y{ ! )[  (ST  i+, )  "  (si-|,i>  1  2 

^(,) - - 5.-1..  h.w'<5.., - - -  ’ 


(2-59) 


where 


S.  .  .  =  f  [  X  ( I )  -  x(i-l)]2  +  [Y(  i)  -  Y(i-l)]2^ 

i  ••  I ,  i 

si,i+i  =  ^x(|+,>  -  +  CYd+i)  -  v(i)]2f 


(2-60) 


The  partial  derivatives  in  the  linearized  stress  displacement  relations, 
Eqs.  (2-58),  must  be  expressed  in  forward  difference  form  at  the  boundary 
of  the  hole  and  are,  for  example, 


u»,  ('.j)  =  ££  O 3 Li ( i ,  j )  +  4u(i+k,j)  -  u(i+2k,j)]  +  0(k2) 
u>2 ( ' » j )  =  C-3u(i,j)  +  4u(i,j+k)  -  u(i,j+2k)]  +  0(k2) 


(2-6l) 


The  linearized  finite  difference  form  of  the  traction  free  boundary  ma^ 
then  be  expressed  as 

iy-  u(i,  j)+X2Au(i+k,  j)  »-;^u(i+2k,  j)+;^Au(i ,  j+k)+;^ju  ( i ,  j+2k) 

+  m?Avl i ,  j  )+m2Av(i+k,j)+m^.‘.v(i+2k,j  )+nyiv(i ,  j+k)+m^iv(i ,  j+2k) 

+  k^  =  0  (2-62a) 


n^u(i  ,j)+n2Au(i+k,j)+ny'u(i+2k,j)+n^u(i,j+k)+n^u(i,j+2k) 

+  P1Av(i,j)+p2iv(i+k,j)+p3iv(i+2k,j)+pZtAv(i,j+k)+p5fi.v(i,j+2k) 

+  k5  =  0  (2-62b) 


where 

l]  “  2k  +  f2*  dS  "  *S1  *  a2*  dS 
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V5k[f>  £  C.J) 


dS 


,  2  dV  .v 

;4  "  k  a2  dS  ^',j) 


S  =  2k  *-f2  dS  *',j) 


a,  g  ('*])] 

f2f  (i.j)] 


d  Y  /  r  •  \  *| 

a2  dS  ("j)1 


■|  *  2k  t(9|  +  %>  §  (i’J>  -  (b1  +  b2>  3? 

”3  "5k[9|  35  (i'j>  -  bi  33 

2  <•,  d Y  >.  dX  i , 

"4  *  T  lbo  jc  ('»J>  "  92as 


"4  k  L  2  dS 

1  r  d> 

n5  =  2k  1 92  dS 


1  r  dX  /  •  •  \  i  d  Y  / .  .  \  ■* 

m-  =  oC  L  9q  ,)<;  (  1  '  J  '  "  b2  dS  1  '  ’  -* 


ki  3?  (i’j>  -  sf  (i'j) 


'1  "  2k 


dX 


dY 


n»  -  [  (dj  +  d2)  ~  ( i  *  j )  -  (f,  +  f2)  ^  (i,j)] 


2  r t  ^  f m  ' \  a  ^  /  *\i 

n2  =  k  [fl  dS  (',j)  "  di  m 


1  dS 


n3  2k  fdl  dS  "  fl  dS 

"4'!tf2  3l  C'J)  -d2S  <'-J» 


"5-k[  d2  3l  (i’j)  -  f2  33  (i-i» 
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P,  '  h  t(el  +  e2>  3?  "  <9I  +  92>  3? 

P2  =  |t9,  f  (>.j)l 

lr  dX  /.  dY  ,. 

p3  =  2k  el  3s  (,'j)  '  9I  d?  (,’J)1 
kt92dS  (,*J>  'e2dS 

1  r  dX  /.  .v  dY  /.  •  \n 

p5  *  2k  ^*2  dS  %  dS 


k5  =  k3  dS  ^  t»s  »-*) 


(2-63) 


j  Y  j  y 

The  values  of  —  (i * j ) »  IT  ( '  » j )  at  the  center  of  each 
db  db 

boundary  cel  1  are  obta i ned  by  averaging  the  values  of  the  derivatives 
of  all  of  the  boundary  particles  lying  in  the  cell  at  any  step. 

Again,  it  should  be  noted  that  the  boundary  conditions  at  each  incremental 
step  are  satisfied  in  the  same  sense  as  they  are  in  classical  elasticity 
theory . 

G.  Extension  of  a  Thin  Elastic  Sheet  with  a  Rigid  Inclusion  at  its  Center 

We  consider  a  rectangular  sheet  subjected  to  biaxial  extensions, 
and  having  a  rigid  inclusion  at  its  center.  The  inclusion  is  of 
arbitrary  shape;  however,  since  we  desire  to  consider  only  one  quarter 
of  the  sheet,  symmetry  of  the  inclusion  is  required.  The  previous 
formulation  for  the  traction  free  hole  remains  unchanged  except  for  the 
boundary  conditions  at  the  hole,  Eqs.  (2-57),  which  is  now  replaced  by 


the  statement  that 
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u  =  v  *  0 

at  the  interface. 

H.  The  Computer  Program  and  the  Computation  Procedure 

The  program  described  here  was  written  in  fortran  language  and 
run  on  an  IBM  360,  model  5 0 >  computer.  In  its  present  form  it 
represents  the  basic  program  necessary  to  solve  the  problems  described. 

No  auxiliary  storage  facilities  are  used  and  the  program  is  presently 
limited,  because  of  the  machine  storage  capacity,  to  solutions  of  small 
finite  deformations.  In  order  to  be  able  to  exceed  the  present  range 
ft  is  necessary  to  include  enough  auxiliary  storage  units  so  that  it 
will  be  possible  to  include  the  new  cells  which  are  required  as  the 
deformation  progresses.  It  is  also  necessary  to  use  a  machine  which 
retains  more  significant  figures  since  the  loss  of  significant  figures 
is  quite  rapid  in  solutions  of  large  matrices  such  as  those  encountered 
here . 

The  program  consists  of  a  main  program  and  a  set  of  subroutines 
The  main  program  reads  in  data,  prints  output,  and  controls  the  sub¬ 
routines.  The  parameters  of  the  problem  are  read  into  the  main  program 
by  data  cards  and  can  be  changed  as  desired.  These  parameters  consist 
of  the  spacing  k,  the  initial  dimensions  of  the  sheet  and  its  material 
coefficients,  the  shape  of  the  hole,  the  size  of  displacement  increments 
and  the  number  of  incremental  steps,  and  the  spacing  of  the  distributed 
particles . 
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The  computation  proceeds  by  first  calling  for  the  subroutine 
which  distributes  the  particles  throughout  the  sheet  and  along  the 
boundary  of  the  hole.  This  step  is  followed  by  a  scanning  subroutine 
which  scans  the  region  of  the  hole  and  determines  the  location  of  the 
boundary  cells,  which  are  defined  as  cells  which  contain  one  or  more  of 
the  boundary  particles.  The  main  program  next  calls  for  the  subroutine 
which  calculates  the  derivatives  ~  ,  — r  at  each  of  the  boundary 
particles  and  then  attributes  to  the  boundary  cell  in  which  they  lie. 
their  average  values.  This  is  then  followed  by  a  subroutine  which 
calculates  the  coefficients  of  the  linear  algebraic  simultaneous  equations 
arising  from  the  finite  difference  solution  of  the  piecewise  linearized 
forms  of  the  equilibrium  equations  (2-49)  and  the  boundary  conditions 
(2-&).  The  same  subroutine  also  calculates  the  stresses  at  the  end 
of  each  step. 

Both  an  iterative  and  direct  method  were  used  to  solve  the 
simultaneous  equations.  The  iterative  method  used  is  the  S.O.R.  which 
can  be  described  briefly  as  follows.  Consider  an  nxn  matrix  A=  (a.j) 
which  is  nonsingular  and  with  nonzero  diagonal  elements.  We  seek  the 
solution  of  the  system  of  linear  equations 


Ax  =  K 


E  a. .  x.  =  k. 

>1  ,J  J  ' 


l  <  i  <  n 


(2-64) 


The  S.O.R.  iteration  process  is  given  by  the  following  algorithm 


*!”"')  .  *!m)  +«(.  -t-  f  a..x.("tl>-  -t-  I 
i  i  a . .  .  .  i  j  j  a.... 

1 1  j=l  J  J  i  i  j=  i 


v  (m)  ■  i  (m). 

E  a. .x.  +  -  x.  \ 

■j-I+1  'J  J  3  i  ?  ‘ 

(2-65) 


where  id  is  the  relaxation  factor  and  (m) ,  (m+1)  indicate  the  m 

and  iteration  respectively.  When  cd  =  1  the  S.O.R.  method  is 

identical  with  the  Gauss-Seidel  iteration  method. 

When  iteration  is  used,  parameters  such  as  maximum  er^or  and 

number  of  iterations  allowed  have  to  be  read  in.  Starting  with 

assumed  initial  distribution  of  Au,  Av,  the  iteration  is  carried  out  by 

moving  along  pivotal  points  which  lie  along  lines  parallel  to  one  of  the 

symmetry  axes.with  Au  calculated  first, fol lowed  by  the  calculations  of 

Av  at  the  same  pivotal  point.  Proceeding  from  the  outer  boundary  towards 

the  hole  or  in  the  opposite  direction  gave  similar  rates  of  convergence. 

The  direct  method  used  is  based  on  the  Gauss  elimination  with 

[39] 

pivotal  condensation  .  First  we  select  the  largest  element  of  the 
coefficient  matrix  A,  say  a^.,  which  is  called  the  pivot.  A  multiplier 
m.  =  -  a. ./a.,  is  then  computed  for  each  row  i^k.  Next,  the  kth  row 

i  i  j  k  j 

is  multiplied  by  m.  and  added  to  the  i^  row.  This  leads  to  a  new 

matrix  with  zeros  in  the  j**1  column,  except  for  a^j  .  The  same  pro- 

th 

cess  is  now  repeated  for  the  auxiliary  matrix  which  excludes  the  j 
column  and  the  k^  row.  The  process  yields  an  upper  diagonal  matrix 
which  can  be  written  with  unit  entries  along  its  main  diagonal.  The  solu¬ 
tion  is  obtained  by  going  from  the  last  row  of  the  upper  diagonal  matrix 


to  the  first  and  in  each  row  substituting  all  the  known  values  and  trans- 
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ferring  all  terms  off  the  main  diagonal  to  the  right  hand  side.  Thus 
the  values  of  the  right  hand  side  yield  the  solution.  The  disadvantage 
of  this  procedure  is  the  rapid  loss  of  significant  figures  which  occurs 
for  large  matrices. 

The  IBM  GELB  subroutine  was  used  to  carry  out  the  direct  method. 
This  subroutine  is  written  for  solution  of  banded  matrices.  It  uses 
Gaussian  elimination  with  col umn  pivoting  only.  The  subroutine  GELB 
requires  that  the  elements  of  the  matrix  be  stored  row  after  row  as  a 
one  dimensional  array.  Therefore,  GELB  is  preceeded  by  a  subroutine 
which  arranges  the  elements  of  the  matrices  in  the  proper  order. 

The  solution  of  the  simultaneous  equations  yields  &u  and  Av 
of  each  cell  center.  Next,  a  subroutine  is  called  for  which  computes 
the  total  displacements  of  each  particle,  their  new  coordinates  and 
cells  in  which  they  lie.  The  new  boundary  cells  are  then  found,  and 
the  magnitude  of  the  displacements  associated  with  the  cell  centers  are 
recalculated.  This  completes  the  computations  for  the  first  increment. 
Further  steps  are  computed  by  reading  in  the  new  prescribed  end  dis¬ 
placement  increments  and  repeating  the  general  procedure  described  above. 

The  program  was  also  adopted  to  the  solution  of  the  problem  of 
the  rigid  inclusion  in  a  sheet.  it  contains  all  the  elements  of  the  hole 
program,  with  the  exception  that  no  derivatives  are  calculated  along  the 
interface  and  that  the  boundary  cel  Is, of  course,  remain  fixed. 
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3.  EXPERIMENTS 


A,  Experimental  Equipment 

A  biaxial  testing  apparatus,  designed  by  the  author,  was  used 
to  conduct  the  biaxial  experiments,  it  consists  of  a  loading  mechanism, 
an  environmental  chamber  and  control  and  measuring  instruments  (Tig.  3). 
The  loading  mechanism  consists  of  an  electric  motor  (1/2  H.P.)  connected 
to  four  shafts  through  a  conical  gear  transmission.  The  shafts  are  thus 
driven  simultaneously  and  have  identical  speeds.  Each  shaft  is  connected 
to  two  parallel  driving  screws  by  means  of  sprockets  and  chains,  resulting 
In  movement  of  the  four  crossheads,  which  in  turn  transmit  their  motion  to 
the  specimen  through  a  loading  frame.  The  stretching  mechanism  thus 
ensures  that  the  center  of  the  specimen  is  always  stationary  and  there¬ 
fore  always  remains  centered.  Crosshead  speeds  in  the  two  principal 
directions  can  be  varied  independently  by  adjusting  the  gear  ratios  in 
the  chain-sprocket  assembly.  A  range  of  crosshead  speeds  from  1  to  20 
in/mfn  can  be  obtained  by  a  suitable  adjustment  of  the  sprocket  assembly 
and  motor  speed.  Each  crosshead  has  a  travelling  range  of  15  in. 

The  test  specimen  (a  rectangular  thin  sheet)  was  cut  so  that 
a  series  of  lugs  was  formed  on  each  of  Its  sides.  Each  of  these  lugs 
was  connected  to  a  light  aluminum  end  strip  which  is  pinned  to  a  roller, 
consisting  of  a  miniature  ball  bearing  mounted  in  an  aluminum  U-section. 
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1  he  rollers  are  constrained  to  rove  along  the  rails  of  the  aluminum 
loading  frames  (Fig. 4).  each  side  of  which  is  connected  to  the  crosshead 
by  twin  parallel  rods  which  pass  through  special  low  friction  guides. 

The  lengths  of  the  rods  are  adjustable  for  use  with  different  size 
specimens.  Out  of  plane  or  lateral  motions  are  prevented  by  support¬ 
ing  the  loading  frames  on  rollers  which  rest  on  an  adjustable  support 
platform.  Loads  are  measured  by  placing  load  cells  between  the  guide 
rods  snd  the  crossheads. 

The  environmental  chamber  Is  an  insulated  box,  having  test 
chamber  dimensions  of  30  *  30  x  12  inches.  Cooling  coils  and  strip 
heaters  are  attached  along  the  side  wails  and  four  sms’?  fans  are  used 
for  air  circulation,  so  as  to  maintain  a  uniform  temperature  throughout  the 
chamber .  The  cooling  is  done  by  connecting  the  coils  to  a  1  ?/2  H.P.  Bendix- 
Westinghouse  Co.  low  temperature  compressor,  and  supplementing  this  by 
using  liquid  nitrogen  for  the  very  low  temperatures.  A  range  of  tem- 
parature  from  -80  to  200°F  can  be  achieved  with  variation  of  ~  1°F 
throughout  the  test  chamber.  Temperature  is  controlled  by  a  West  instrument 
Corp.  Gardsman  control  unit  which  controls  the  strip  heaters.  Below 
room  temperatures  are  reached  by  operating  the  compressor  continuously 
and  controlling  the  heat  input  of  the  heaters  It  was  found  that  this 
operating  procedure  results  in  a  more  stable  temperature  than  when  both 
heating  and  cooling  rates  or  only  cooling  rates  are  controlled.  The 
upper  surface  of  the  environmental  chamber  is  covered  by  e  thermal  pane, 
free  to  move  along  two  support  rails,  thereby  permitting  ootical  measure¬ 
ments  of  the  displacements  of  the  test  piece,  as  well  as  easy  access  to 


the  chamber  interior. 


Sane  uniaxial  experiments  were  also  conducted.  -  The- equipment 
used  for  these  tests  consists  of  a  Universal  Testing  Instrument  Model 
TTC  Hi .3  made  by  Instron  Engineering  Corporation,  and  an  FRL  Environ¬ 
mental  Box  designed  by  Fabric  Research  Laboratories.  Heating  in  the 
chamber  is  accomplished  by  use  of  strip  heaters,  and  cooling  is  done 
by  use  of  either  dry  ice  or  liquid  nitrogen.  A  Barber-Colman  *'o . 
Wheelco  unit  provides  temperature  control  over  a  range  of  -100  to  400°F 
by  operating  the  strip  heaters. 

8.  Specimen  Preparation 

The  experiments  were  performed  on  specimens  of  natural  rubber 
(cured  with  1  part  of  di  cumy!  peroxide,  2.5  pt  dicup  40c)  all  of  which 
were  cut  from  the  same  sheet  in  order  to  reduce  possible  variations  of 
the  material  properties.  The  thickness  of  the  sheet  was  0.079  ±  0.004 
in,  A  specimen  having  6  1/2  in.  sides  was  used  for  the  biaxial 
experiments,  and  a  strip  having  a  length  of  5  in.  and  a  width  of  0.375 
in,  was  used  for  the  uniaxial  tests. 

For  the  biaxial  experiments  a  silk  screen  was  prepared  having 
sub-divisions  varying  from  1/4  in.  at  its  perimeter  to  1/32  in.  in  its 
interior.  Although  the  small  sub-divisions  were  not  required  for  the 
characterization  experiments,  they  were  necessary  for  the  stress  concen¬ 
tration  experiments  which  were  carried  out  on  the  same  sheet. 

The  grid  was  printed  on  the  specimen  b/  placing  the  screen 
over  it  and  spreading  a  rubber  base  ink  (Ridgway  RC  1154,  R.l.  White) 
over  the  surface  of  the  screen  by  pressing  a  squeegee  over  it. 
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C .  Experimental  Procedure 

The  specimen  was  placed  in  the  loading  rig  and  stretched  to  a 
predetermined  elongation  in  one  direction.  This  elongation  was  then 
kept  fixed  while  the  perpendicular  sides  were  stretched  over  a  range  of 
elongation  ratios  from  1.0  to  1.8  in  0.1  steps.  The  procedure  was 

repeated  for  different  values  of  the  fixed  elongation.  Isotropy  in  the 
undeformed  state  was  checked  by  repeating  the  tests  with  interchanged 
principal  axes,  and  noting  that  there  was  only  a  slight  variation  in 
loads  corresponding  to  interchanged  elongations.  Before  each  test  the 
specimen  was  relaxed  for  two  hours  et  the  highest  temperature  of  the 
tests  (150°F)  in  order  to  release  possible  residual  stresses. 

The  experiments  were  performed  at  30,  60,  SO,  120  and  150°F. 

Only  slight  temperature  effects  could  be  detected  in  the  biaxial 
experiments,  and  since  these  were  within  the  experiment.?!  e-rrer, 
temperature  effects  could  not  be  isolated. 

The  displacements  of  the  sheet  were  recorded  s*  prescribed 
cross head  positions  by  photographing  the  specimen  with  a  4  X  In.  view 
camera  (made  by  Calumet  Manufacturing  Co.)  mounted  on  a  tripod  placed 
on  the  thermal  pane  cover  (Fig.  5)-  By  using  high  contrast  film 
the  image  on  the  negative  appears  os  sharp  black  lines  on  a  clear 
background.  A  typical  negative  is  shown  ir  Fig,  6.  The  scales  of 
the  nogativec  were  determined  by  placing  a  standard  ruler  in  the  plane 
of  the  specimen  while  photographing .  Distances  between  grid  lines 
were  determined  by  scanning  the  negatives  with  a  microdens i tome ter 
(made  by  Joyce,  Loebi  and  Co..  England)  or  an  x-y  micrometer  table  mounted 
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on  a  precision  transmission  poia.' scope  (made  by  Photoelastic  Inc.). 
Measuring  sensitivities  of  both  instruments  were  at  least  of  the  order 
of  .0001  in.  The  undfcformed  specimen  was  similarly  analyzed  and  thus 
displacements  and  elongation  ratios  were  determined. 

The  biaxial  tests  were  supplemented  by  a  series  of  uniaxial 
tests.  A  strip  cut  from  the  same  sheet  was  tested  st  temperatures 
varying  from  30  to  150°F.  Temperature  effects  in  the  uniaxial  tests 
were  more  pronounced  than  those  observed  in  the  biaxial  experiments. 

After  completion  of  the  material  characterization  experiments 
a  0.5  in.  diameter  hole  was  punched  in  the  center  of  the  sheet.  A 
series  of  st""-.*:  OTfrencrat Ion  experiments  were  conducted  at  90°F» 
using  same  exte;*  >vsnt.al  procedure  as  described.  A  typical  photo¬ 
graph  of  the  deformed  specimen  is  shown  in  Fig.  7. 

0.  Material  Characterization  experiments 

i)  Analytical  Analysis 

Consider  a  thin  rectangular  plane  sheet  of  uniform  thickness 
having  edges  parallel  to  the  x  and  •/  axes.  The  sheet  is  subjected 
to  a  pure  homogeneous  deformation  having  principal  directions  parallel 
to  the  x,  y  and  z  axes 

Let  ,  ^2»  ^  the  principal  extention  ratios,  defined  as 
the  ratios,  in  the  principal  directions,  of  final  length  to  original 
lengths  of  the  elements 
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The  sheev  is  assumed  to  be  isotropic  and  incompressible,  so  that  the  con¬ 
stitutive  equations  can  be  expressed  as  follows^  ^ 


r  -  2  ('  ^  ^ 
C  5  5  2  (M 


1  DW  . 

p 


(3-2) 


(i  =  1,2,3  or  x,y,z) 


where  c..  is  the  principal  true  stress  and  p  is  an  arbitrary  hydro- 


i  i 


static  pressure.  The  strain  energy  function  W  is  expressible  as  a 
function  of  the  invariants  Tj,  7 which,  in  terms  of  the  principal 
extension  ratios,  are 


p  p  p  p 

Tl  *  V  Kl  *  UX?'2 


'  2  '  +  IAf  +  T  '2 


(3-3) 


Since  c ■  =  0  we  can  solve  for  p  and  rewrite  Eqs.  (3-2)  as 


c22  =  2('|-  IA^)f-+^^-) 

I  £_ 


Solving  for  §7—,  from  Eq .  (3-4)  yields 
*  1  aT2 


(3-4) 
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ii)  Experimental  Results 

The  biaxial  tests  were  conducted  by  keeping  one  of  the  extension 

ratios  fixed  and  varying  the  other.  From  these  experiments  the  stresses 

were  determined  as  functions^of  the  extension  ratios,  and  »  and 

thus  as  functions  of  the  invariants,  Tj  and  T^.  These  values  were 

ay  aw 

used  in  Eqs .  (3-*>a)  and  (3-7)  to  determine  rr—  and  5=—  as  functions 

0- 1  QT2 

aw  1  ^ 

of  T,  and  T_.  Plots  of  r=—  as  a  function  of  T  and  are 

Id  0I1  Id 

shown  in  Figs.  8.  It  is  seen  that  within  the  experimental  error 

*T1 

is  independent  of  both  Tj  and  T^,  and  is  equal  to  the  constant 
Cj  =  20.28  psi.  Similar  plots  of  as  a  function  of  Tj  and 

are  given  in  Figs.  9  *  In  Fig.  9a  we  observe  that  a  reasonable  fit 

aw 

to  the  data  is  that  is  independent  of  1^,  and  Fig.  9b  shows 

that  its  value  decreases  with  increasing  values  of  1^.  A  least 
square  parabolic  fit  to  the  data  yields  the  function: 


5Y-=  [h.808  -  1.440  (I2  -  3)  +  0.1379  (I2  -  3)2]  psi 


(3-10) 


which  is  shown  as  a  solid  curve  in  Fig.  9b.  Thus  we  may  express 
the  strain  energy  function  as 


W  =  20.28(1  j  -  3)  +  5.808(I2-3)  -  (r~-3)2  +  ^12  (T2-3)3  (3-11) 


It  should  be  noted  that  this  form  has  been  obtained  from  data  in  the 
low  ranges  of  Ij  and  I?,  i.e.,  1^  and  I2  <  8. 

Since  the  stress  concentration  experiments  will  be  conducted 
so  that  the  maximum  values  of  and  I2  are  within  this  range,  these 

low  range  results  will  be  applicable. 


Uniaxial  extens <on  tests  were  conducted  at  30,  6^,  90,  120 
and  150°F.  The  results  show  (Fig.  10)  that  the  material  becomes 
stlffer  with  increasing  temperature,  a  result  which  is  in  agreement 
with  the  kinetic  theory  for  rubber  elasticity. 

Fig.  11  shows  a  plot  of  the  uniaxial  stress  as  a  function  of 
elongation  ratio  at  a  temperature  of  90°F.  The  solid  line  has  been 
determined  by  using  the  form  of  the  strain  energy  function  obtained 
from  the  biaxial  experiments  and  substituting  it  into  Eq.  (3-9). 


4.  NUMERICAL  CALCULATIONS 


The  modified  P.I.C.  method  was  applied  to  the  problems  of  a 
circular  hole  and  a  rigid  inclusion  In  a  sheet.  Preliminary  calculations 
were  performed  on  the  former  problem  by  assuming  a  4  x  4  unit  square 
sheet  having  a  circular  hole  of  radius  =  1  unit.  A  2o  x  2o  grid 
(k  -  0.20)  was  used,  thus  resulting  In  approximately  800  simultaneous 
equations  (u,v  are  unknown  at  each  cell  center).  The  equations  were  so 
arranged  that  a  diagonally  banded  matrix  containing  over  100  elements  was 
obtained.  Both  direct  and  iterative  methods  were  applied.  Difficulties 
with  convergence  of  the  S.O.R.  iterative  scheme  arose,  while  In  the  direct 
method  (Gaussion  elimination  with  column  pivoting)  a  great  many  significant 
figures  were  lost  after  several  incremental  steps.  While  results  showed 
agreement  with  the  linear  theory  for  the  first  few  incremental  steps, 
no  meaningful  results  could  be  obtained  for  finite  deformations. 

For  the  rigid  Inclusion  problem  the  displacement  boundary  conditions 
did  not  result  in  a  complication  of  the  matrix,  as  did  the  stress  free 
boundary  conditions  of  the  former  case.  A  10  x  10  grid  (k  =  0.3)  was  used, 
and  the  direct  method  of  solution  yielded  reasonable  results.  Since  the 
displacement  functions  were  much  smoother  for  the  case  of  the  rigid  inclusion 
than  for  the  case  of  the  cutout,  it  is  to  be  expected  that  a  courser  grid 
could  be  successfully  applied  to  the  former  problem. 

It  should  be  pointed  out  that  the  analytical  development  of  the 
linearized  equations  of  motion  (2-47)  was  carried  out  under  the  assumption 


that  the  material  could  be  described  by  the  Rivl in-Mooney  constitutive 
equations.  However,  the  results  of  the  experiments  on  the  natural  rubber 
indicate  that  the  materia!  is  described  by  Eq„  (3—1 1 ) *  and  thus  will  be 
referred  to  as  a  General ized  Rivl in-Mooney  material.  To  make  these  two 
results  compatible,  it  was  decided  that  in  the  analytical  work,  would 
be  chosen  to  be  a  function  of  as  shown  in  Fig.  (9b),  that  is,  it  was 
assumed  that  =  <5W  (I^) /ST^ 


5.  DISCUSSION 


The  results  of  the  material  characterization  experiments. on  a 
thin  sheet  of  natural  rubber  are  presented  in  Figs.  (8-1 1 ).  Figs . (8)  show 
that  Is  independent  of  the  strain  invariants  1^,  I^,  and  has  a  value 

Cj  =  20.28  psi.  Similar  graphs  for  are  found  In  Figs. (9)*  which 

indicate  that  is  independent  of  and  decreases  with  increasing 

values  of  I,,.  It  should  be  noted  that  the  form  of  the  strain  energy 
function  (3-11)  has  been  obtained  from  data  in  the  low  ranges  of  and 
I ^ i  i.e.  for  Ij,  <  8.  Since  the  stress  concentration  experiments  were 
conducted  so  that  the  maximum  values  of  Ij  and  fall  within  this  range, 
these  low  range  results  are  applicable. 

At  the  completion  of  the  material  characterization  experiments 
a  circular  hole  was  cut  at  the  center  of  the  sheet.  A  series  of  biaxial 
experiments  were  then  performed  on  the  specimen  3t  a  temperature  of  90°F. 
Fig.  12  shows  the  experimentally  obtained  displacements  along  the  axis  of 
symmetry  (x=o)  for  the  equi-biaxlal  elongations.  Investigation  of  the 
slopes  In  the  vicinity  of  y/r^  =  13.0  Indicates  that,  essentially,  an 
isotropic  strain  state  exists  there.  Thus  the  results  obtained  ere  not 
very  different  from  those  that  would  be  obtained  from  an  infinite  sheet. 

It  Is  for  this  reason  that  the  ratio  of  the  deformed  to  the  undeformed 
length  of  the  membrane  is  designated  \w(  =  X^/x^  *  Yj/yj)  for  the 
equi-biaxlal  case.  For  the  general  biaxial  state  (\X)M  and  (\  )w  refer  to 
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the  length  ratios  in  the  x-  and  y  directions,  respectively.  Displacements 

for  the  uni-'d irect ional  extension,  i.e.  C'  )  =1.0,  are  presented  in 

x  ® 

Figs.  13,  and  the  corresponding  deFormeu  shapes  of  the  hole  are  shown  in 
Fig.  14.  Substituting  the  displacements  and  their  derivatives  into  the 
constitutive  equations  (2-36)  results  in  the  experimentally  determined 
stresses.  The  stress  ~YY  (circumferential)  at  the  intersection  of  the 
hole  and  the  x  -  axis  is  plotted  in  Fig.  15.  A  comparison  with  linear 
theory  indicates  that  no  significant  deviation  occurs  untii  *.  >1.3. 

A  plot  of  the  stress  concentration  factor  Is  presented  in  Fig.  16, 
where  it  is  seen  that  the  factor  varies  from  ?.Q  for  small  strains  to 
approximately  4.0  at  V  =  1.6. 

No  analytical  results  were  obtained  for  the  problem  of  the  hole 
in  the  sheet  for  reasons  which  are  discussed  si>  Chapter  IV.  Although  no 
experiments  were  performed  on  the  sheet  containing  a  rigid  inclusion,  the 
modified  Particle  -  in  -  Cell  (P.f.C.)  method  was  applied  to  the  problem. 
The  results  for  the  equi-oiaxial  case  are  presented  in  Figs.  17-20.  The 
stress  concentration  factor  K^,  is  presented  in  Fig.  20  for  a  Generalized 
Rivl in-Mooney  (3-11),  a  R  i  vl  i  n-Mooriev  (2-33)  and  a  linear  elastic  material. 
The  curves  indicate  that  the  stress  concentration  factor  decreases  with 
an  increase  in  \  for  the  Generalized  R ivl in-Mooney  material,  whiie  for 
a  Rivl in-Mooney  material  the  situation  is  reversed.  The  results  for  the 
R ivl in-Mooney  material  agree  qualitatively  with  those  obtained  by  Yang  [32], 
and  with  those  obtained  from  an  analysis  of  the  displacements  presented  by 
Rivl In  and  Thomas  [31]. 
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APPENDIX  A 

COEFFICIENTS  OF  LINEARIZED  EQUILIBRIUM  AND 
STRESS-DISPLACEMENT  EQUATIONS 

Let  the  displacements  at  the  end  of  the  n-th  step  be  denoted  by  u 
and  v.  The  coefficients  of  the  linearized  equations  at  the  (n  +  l)-th  step 
are  then  given  by  the  following  expressions, 

i)  Coefficients  of  Linearized  Equilibrium  Equations  (2-47) 

Aj  =  -C  [2  V(1  -u^)  +  2  U  u,2  +  3  L4  v>12  +  12  L3  M  (1  -  v>2)] 

+  L  {-  2L  [3(1  -u.jHl  -v,2)  *  u,2  v. -j  ]  u,j,  -  2S  (i  -v,?)  u.jj 

-I"  2L  [(1  -  v»  2)  u»2  "  ^  j  j  \2.~  "^>2)  u*  12“  2L  ^ ^ 

-2Z  (1  -v,2)  u,22  +  2L[{1  -v,2)  V>1  -  (I  -u,j)  u,2]  v,n 
■  2 P  (1  -v,2)  v,n  +L[(u>2)2-  3  (1  -u,j)2-  (v,j)2-  (1  “V,2)2]  v,12 

-  2Q  (1  -  v,2)  v,12+  2  LT  v,22  -  2R  (1  -v>2)  v>22  }  -  v>12 

A2  =  C[Wv,22+2Tv,12  +  2Y  u,2  +  2U  (1  -u.j)  +  XV>11  +  3L4  v,n-12L3Mvfl] 

+  L  {-  2L  (1  -u.^tu.jj  +  u,22)v,1  -  2L2u,12  -  2Nv,,  u,12- 2Su>n  v,j  -2  Z  v^u^j 
+  L  [W  +  2(v,  jl2  ]  v,  x  j  -  2Pv,  jV,  jj-2  LT  v>12~  2Qv>12  v,j  +  L  [(v^)2 

-  (!  -  u.j)  ]  v,22  -2R  v,j  v,22  }  +  v.jj 

Aj  j  =  C  [  X  (1  -  v,2)  +  3  L4  (1  -  v,2)  ]  +  L2  S  +  (1  -  v,2) 

A12  =  C  [  2T  (1  -  v,2)  +  3  L4  v,j]  +L2  N  +  v,j 
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Az?  =  C  W  (1  -  v,2)  +  L2  Z 

B,  =  C  [2  V  v,j  +  2  U  (1  -  v,2)  +  3  L4  u,12  -  12  L3  M  u,2  ] 

+  L  {  2L  [(1  -  v,2)  v,j  -  (1  -u.j)  u, 2  ]  u,n  -  2  S  u,n  u>2 

+  L  (  (1  -u.j)2  +  (u,2)2  +  3  (v,j)2  -  (1  -v,2r  ]  u,12  ‘  2LT  u>22  -2Zu,22u>2 

+  2L[  3  v,  j  u,  2~  (1  -utl)(l  -v,2)]  v,jj-  2  Pv.jj  u,2  -  2Nu,J2  u>2 

2 

Hf-  2L[(1  * u.,  j)  V|  j-  (i  •  vf  ^  1 2 *  2  ^  12*  ^  22 *  ^  2  22  ^ 

+  u'12 

B2  =  -  C  [Wu,22+  2Tu,12+Xu,n  +  2Y(l -v,2)  +  3L4u,n+2Uv,1+12L3M(l-u,1)] 

+L  {  -L  [  3  (1  -u,1}2  +  {v,1)2]u,I1-2S(l  -u11)u,11+2L,Tu,12-2N(1  -u,j)u,12 
+  L  [(v.j)2  -(1  -utl)2]u,22-  2Z{1  -u>;l)u,22+  2L  (1  (v»11+v»22) 

-  2P  (1  -u.j)  v>u  -  2L2  v>12  -  2Q  (1  -  u,j)  v,12 

-  2R  (1  -  u,j)  v>22  }  -  u,  u 

Bn  =  C  [X  u>2  +  3  L4  u,2]  +  L2  P  +  u,2 

B12  =  C  [2T  u,  2  +  3L4  (1  -  u>1)  ]  +  L2Q  +  (1  -  u,j) 

B22  =  CWu,2  +  L2P. 

Kj  =  C  [WV  +  2TU  +  XY  +  ?  L4  M  ]  (A-l) 

+  Lu  [S  u.jj  +  N  u,  12  +  Z  u,22  +  P  v.jj  +  Q  v,J2  +  R  v>22]  +  M 


where 


C  -  Cj/C2 


L  = 

(1  - 

■  u,  j)(i  -  y>2)  -  u>2  v,j 

M  = 

0- 

v»2)  u»  +  (1  -u,  j)  v,  12 

+  11,21  v, 

1  +u-2 

v’ll 

N  = 

(1  - 

u,j)2  y.j  +(v,1)3-2u,2 

(1  -ufl)(l 

-v,2)- 

2  2 

v,  i  (!  ~v,2)  +(u,2)  vfl 

P  = 

a,2 

(1  -  u,  j)2  -  2  vfl  (1  -  u. 

1)(1  "  v,2 

)  +  3  u. 

?.  (vM)2 

Q  -• 

(1  - 

u,j)3  +  (1  -u.jlKv.j)2  - 

(u*2)21  " 

2u>2  v. 

j  (1  -v,2)  +  (l  -v,2)2(1  - u, j ) 

R  = 

U’2 

(v,j)2  -  u>2  (1  -  u,j)2  - 

2v,j  (1  - 

■  u,  j)(l 

-  v,2) 

S  =  (1  -  v^Kv.j)2  +  3  (1  -  u^)2  (1  -  v>2)  -  2  (1  -  u,j)  u^v.j 

T  =  u,2  (1  -  u,j)  +  v,j  (1  -  v,2) 

U  =  (!  ~  v,2)  ix,i2  +  n,2  v»i2 

V  =  (1  -v,2)u,22+u,2  v,22 

W  *  (1  -  u.j)2  +  (v.j)2 

X  =  (1  -  v,2)2  +  (u,/ 

Y  =  (1  '  V'2)U>11  +  u’2  V’ll 

2  2 

Z  —  (1  —  U|  j )  (1  "  v>2)  _  2  u>2  |  (1  ~ u,  j )  ~  \v ,  j)  (1  —  ^’2'  (A  ~2) 

and 

Dj  =  -C  [Wv,n  +  2T  v,12  +  2Uu?2  +  Xv,22  +  3E4v,22  +12  L3M(1  -v,,) 

+  2  Y  (1  -u,j)]  +L{2L(1-v,2)u,2(u,11  +  u,22)-2R  (1-v,2)u,  jj  -  2L2  u,,2 
-  2Q(1  -  v,2)  u,12  -  2P  (1  -v>2)  u>22  +  l.[(u.,2)2  -  (1  -  v,2)2},v,  n 


>4?. 

-  2Z  (i  -v,2)  V.JJ  +  2LT  v,12  -  2N  (1  -v,2)  v,  12~  L,[(u,2)“+  3{1  -v,  2)Z]  v, 
-2S(1  -  v,2)v,22  }  -  v,22 

D2  =  C[2Vu)2T2U(l-u,1)+3L4v,12-12L3MvI1] 

+  L{=*2L2u,  j  j  -  2  R  v,  j  u.jj  +  2  L[(l  -  v>2)  u>2  -  (1  -  u,j)  v,j]  u,12 

-  2Qv-1  u,12  +  2  L  [3  u,2  v.j  -  (1  -  u,j}(l  -v>2)]  u,22  -  2  Pu,^  v,j 

-2  LT  v.jj-2  Z  vtll  v,,  +L[(1  -v,2)2  +  3  (u>2)2  +  (v,  j)2- (1  -u,j)2]  v,12 

■  2  N  v,  j  v,  ^  +  2L  [u,  2  (1  _u»  j)“  v>  j  ~v,2^v,22  ”  ^  ^  V,ZZ  v*  1  ^ 

+  v,12 

Dn  =  C  W  V>1  +  L2  R 

D12  =  C[2T  v,j  +3L4(1  ■  v, 2)]  +  L2  Q  +  1  -  v>2 
d22  =  C  [  X  vfl  +  3  L4  v.j]  +  L2  P  +  V>1 

E1  =  C  [  W  u,n  +  2  Tu,12  +  X  u,22  +  2  YvfJ  +  3  L4  u,22  -  12  L3  Mu,2 

+  2  0  (1  -  v>2)J 

—  2  2  ••  _ 

+  L  {-  L[{1  -  v,2)  -  (u>2)  ]  ufJ1  -  2R  ufll  u,2  -  2LT  u,12 

"  2  Q  U'12  u,2  +  L  K1  "  v'2^  +  3<u*2^2  1  u,22  ’  2  P  u»22  u»2 

-  2  L  (1  -  v>2)  u, 2  (v>n  +  v,22)  -  2  Zv,n  u,2  -  2  L2  v,^-  2N  v>l2  u>2 

•2^  v>  22  u'2  }+u’22 
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E?  =  C  [-2V  (1  -  v,2)  -  2  U  vf  j  -  3  L4  u,  12  -  12  L3  M(1  -  u^)] 

+  L{2LT  u,  j  j  -  2  E(1  “  u,  j)  u*  X1  +  L  [(v»  j)2-^,  2)2  -3(1  -v,  2)2 
-(1-u,  j)2]  u,  12  -  28  (1  -u,  )  u,  12+2L[(1-u,  j)  u,  2  -  (1-v,  2 )  v  .]  u,  22 

-  2P  (1-u,  j)  u,22-2L2  v,  n  -  2Z(1  -u,  x)  (v,  n)  +  2“  [(1-u,  j)  v, , 

-  (1  -v,  2 )u,  2]v,  i2  -2N  (1  -u,  j)v,  I2+  2L[u,2  v,  j-3  (1  -u,  j)(l  -  v.2)]v»22 

-  2S  (1  -u,  j)  v,  22>  —  Uj  12 

En  =  C  W  (1  -  u,  +  L2  Z 

E12  =  C  [2T  (1  -  u,  x)  +  3L4  u,  2]  +  L2  N  +  u,  2 

E22  =  C  [X  (1  -u,  L )  +  3L4  (1  u,  j)]  +  L2  S  +  1  -  u,j 
K2  =  C  [  W  V  +  2  T  U  +  XY  +  3L4M] 

+  L2[Ru,11  +  Qu,12+Pu,22+Z  v,  Nv,  12+  S  v>22]  +M  (A-3) 

where 

L  =  (1-u,  j)(l-v,  2)  -  u,  2  v,  l 

M  =  (1-u,  j)  v,  22+  (1  -v,  2)  u,  12+  v,  |2u,  2+  v,  1  u,  22 

N  =  (1-v,  2)2  u,  2+  (u,  2)3  -  2  (1-u,  j)  (1-v,  2)  v,  j-  (1-u,  j)2u,  ^(v.^u,  2 

P  -  V,  J  (1-v,  2)2  -  2  u,  2  (1-u,  j)  (1-v,  2)  +  3  (u,  2)2  V,  J 

8  =  (1-v,  2)3+  (1-v,  2)  [(u,  2)2-  (v,  j)2]  -  2u,  2v,  l  (1-u,  x)  +  (1-u,  ^(l-v,  2) 

R  =  [(u,  2)2-  (1-v,  2)2]  v,  j-2  u,  2  (1-u,  j)  (1-v,  2) 

S  =  (1-u,  j)  (u,  2)^  +  3  (1-v.  7V‘  (1-u,  j)  -  2  (1-v,  2)  u,  2v,  j 
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T  =  v,  1  (1-v,  2)  +  u,  2  (1-u,  j) 

U  =  (l-u.j)  v,  12  tv.ju,  ,2 

V  =  (1-u,  j)  v,  n  +  v,  1  u,  X1 

W=  (1-v,  2)2  +  (u,  2)2 

X  =  (l-u.j)2  +  (V>1)2 

Y  =  (1-u.j)  v,22  +vtlu,22 

Z  =  (1-v,  2)2  (1-u,  j)  -  2u,  2  v,  ^1-v,  2)  -  (u,  2)2  (1-u,  x)  (A -4) 

ii)  Coefficients  of  Linearized  Stress -Displacement  Equations  (2-58). 
ax  =  4{2C1L3(1-v.  2)  +  C2L  [W(1-v,  2)  +  L(l-u,  j)]} 

a2  =  4[Cj(2v,  XL3  +  u,  2)  +  C2L  W  v,  j] 

bj  =  4[2Cj  L3u,  +  C2L  (W  u,  2-  v,  jL)] 

b2  =  4{C1[2L3  (1-u,  j)  -  (1-v,  2)]  +  C2L  W  (1-u, 

kx  =  2[C1(X-L4)  +  C2(l-L2  W)] 

dj  =  4{C1[2L3(1-v,  2)  -  (1-u,  x)]  +  C2L  X  (1-v,  2)} 

d2  =  4[2C1L3  v,  1  +  C2L  (Xv,  1  -  Lu,  2)] 

ex  =  4[C1(2L3u,2  +vfl)  +  C2L  Xu,  J 

e2  =  4{2CX  L3(1-u,  x)  +  C2L  [X  (1-u,  j)  +  (1-v,  2)  L]} 

k2  =  2[CX  (W-L4)  +  C2(l  -  L2  X)] 

fj  =  -2{CX  u,  2  +  C2L  [L  u,  2+  2  T  (1-v,  2)]} 
f2  =  2(^(1-^  x)  +  C2  L  [L(l-u,  x)-2  T  v,  J  } 

gj  =  2{C1(l-v,  2)  +  C2  L  [L(l-v,  2)  -  2  T  u,  2]} 

g2  =  -2{C1  v,  j  +  C2  L  [L  v,  1  +  2  T  (1-u,  j)]} 

k3  -  2  T  [Cj  +  C2  L2] 


(A -5) 


Symmetry  Axis 


FIG.  1  EULER  IAN  GRID 
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|j.  7  PHOTOGRAPH  0"  DEFORMED  GRID  (GRID  SPACING:  1/3?  -  1  /*+  in.') 


FIG.  20  STRESS  CONCENTRATION  FACTOR  AT  RIGID  INCLUSION  (K  ) 

r 

EQUI-B IAX IAL  EXTENSION 


UNCLASS iFI ED 


Security  Classification 


DOCUMENT  CONTROL  DATA  -RID 


_ _ fSKurifr  Ol maelfimetien  ol  title.  body  ol  ebetrect  end  indexing  mnm.jtenon  must  bo  ontoro-f  when  tho  oworotl  report  i»  cloaotfiod) 


\*  ORIGINATING  ACTIVITY  (Corporate  euihor)  I  2a.  REPORT  SECURITY  CLASSIFICATION 

Polytechnic  Institute  of  Crooklyn  |  UNCLASSIFIED 

Dept,  of  Aerospace  Engineering  6-  Applied  Mechanics 
333  Jay  5t.»  Brooklyn,  K.  Y.  11201 


s.  report  title 

STRESS  CONCENTRATION  IN  AN  ELASTOMERIC  SHEET  SUBJECT 
TO  URGE  DEFORMATIONS 


4.  DESCRIPTIVE  NOTES  (Type  of  report  end  inclooioe  dot  mm) 

iNTERIM  TECHNICAL  REPORTS 


5.  AUTHORISE  (Pixel  nemo.  middle  inltiml t  toot  nemo) 


Alexander  Segal  and  J.  M.  Klosner 


C.  REPORT  DATE 


MARCH  1970 


•a.  contract  or  grant  no.  DAHC04-69-C-001 9 


OTAL  NO.  OF  PAGES  176.  NO.  OF  REFS 


9a.  ORIGINATOR’S  REPORT  NUMUERISi 


6.  PROJECT  NO. 


‘  Pibai  Report  #70-11 


*  9b.  OTHER  REPORT  NOlS)  (Arty  otf.er  numbers  the;  nay  be  assigned 
j  thin  report) 


10.  DISTRIBUTION  STATEMENT 


U.  S.  Army  Research  Office 

Distribution  of  this  document  Is  unlimited 


h.  supplementary  notes 


12.  sponsoring  military  activity 


U.S.  Army  Research  Off ice-Durham 


m 


ABSTRACT 

Biaxial  and  uniaxial  experiments  have  been  conducted  on  a  thin  sheet  of  natural 
rubber,  which  can  be  assumed  to  be  incompressible,  isotropic,  and  perfectly  elastic. 
The  strain  energy  function  and  consjitutive  equations  have  been  determined,  and  the 
material  Is  classified  as  a  Generalized  Rivl in-Mooney  type. 

Biaxial  experiments  were  then  conducted  on  the  same  sheet  with  a  circular  cutout 
and  stress  concentration  factors  were  obtained.  Results  indicate  a  significant 
increase  in  the  factor  with  Increased  displacements. 

A  modified  Part icle- ln-Cel  1  (P.jT.C.)  method  has  been  developed  and  analytical 
results  were  obtained  for  a  sheet  with  a  rigid  circular  inclusion.  It  is  shown  that 
the  stress  concentration  factor  for  a  Rivlin-Moor.ey  material  increases  with  increasing 
deformations,  a  result  which  is  in  qualitative  agreement  with  solutions  obtained  by 
other  methods  [31»3S?dK  The  use  of  the  Generalized  R I  vl  i  n-Mooney  material,  however, 
leads  to  a  decrease  in  stress  concentration  with  increasing  deformations. 
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